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ABSTRACT 

Using the SDSS spectroscopic sample, we estimate the dark matter halo spin parameter 
A for <~ 53, 000 disk galaxies for which MOPED star formation histories are available. 
We investigate the relationship between spin and total stellar mass, star formation 
history, and environment. First, we find a clear anti-correlation between stellar mass 
and spin, with low mass galaxies generally having high dark matter spins. Second, 
galaxies which have formed more than ~ 5% of their stars in the last 0.2 Gyr have 
more broadly distributed and typically higher spins (including a significant fraction 
with A > 0.1) than galaxies which formed a large fraction of their stars more than 10 
Gyr ago. Finally, we find little or no correlation between the value of spin of the dark 
halo and environment as determined both by proximity to a new cluster catalogue 
and a marked correlation study. This agrees well with the predictions from linear 
hierarchical torquing theory and numerical simulations. 
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1 INTRODUCTION 

The influence of tidal torques on the evolution of galaxies 
was recognized very early on (Hoyle 1949; Doroshkevich 
1970) as a driver of galaxy morphology. Later on, the spin 
parameter A of dark matter halos, which was introduced by 
Peebles (1969) as 



J\E\ 



1/2 



(1) 



where J is the angular momentum, E is the energy, and M 
is the total mass of the galaxy, was studied analytically and 
numerically within the framework of hierarchical galaxy for- 
mation by several authors (e.g. Barnes & Efstathiou (1987); 
Heavens & Peacock (1988); Warren et al. (1992); Catelan & 
Theuns (1996)). 

Vitvitska et al. (2002), Mailer et al. (2002) and Hetz- 
necker & Burkert (2006) have put forth that in addition to 
tidal torquing, galaxy mergers play a driving role in the evo- 
lution of galactic spins. While they find that recent major 
mergers increase a halo's spin, D'Onghia & Navarro (2007) 
emphasize that mergers influence spin only as long as the 
halo is out of equilibrium, and that isolated, virialized haloes 
have spins consistent with simple tidal torques. 

It is reasonable to expect that the spin of a dark matter 
halo will influence the final properties of baryonic matter in 
the galaxy. For a simple example, consider a disk galaxy and 
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assume baryons settle into the disk with no loss of angular 
momentum. As given in Eq. 1, A measures the degree to 
which rotation contributes to supporting the galaxy against 
collapse, between negligibly (A = 0) and completely (A ~ 
1). Higher A disks are more rotationally supported and will 
therefore be less dense. When coupled with a star formation 
law dependent on density (Kennicutt 1998), higher A further 
implies less efficient star forming systems. 

Using analytical models, Dalcanton et al. (1997), 
Jimenez et al. (1997), van den Bosch (1998), and Mo et al. 
(1998) showed in detail how a distribution in the values of 
the halo spin parameter could lead to significant differences 
in the star formation efficiency of the disk, thus shaping the 
history of the galaxy, and therefore the morphological type, 
even leading to possible cases of dark galaxies (Dalcanton et 
al. 1997; Jimenez et al. 1997; Verde et al. 2002), where star 
formation in the disk has been completely prevented. 

It is of great interest to measure observationally the 
dark halo spin and determine its link star formation, espe- 
cially in view of the possible role of mergers on influencing A 
and their known link to star formation. In a previous study, 
Jimenez et al. (2003) used a sample of observed rotation 
curves of spiral galaxies to determine the spin distribution 
P(A), showing that it was consistent with that predicted by 
hierarchical tidal torque theory (see their Fig. 2). Jimenez 
et al. (2003) also showed that the spin parameter correlated 
with the baryon fraction in the galaxy: the higher the spin 
the higher the baryon fraction. 
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Studies of the spin parameter in dark matter numeri- 
cal simulations have recently reached maturity in that the 
number of particles per halo is now large enough for statis- 
tical samples of the spin parameter to be reliably measured 
(Nagashima & Gouda 1998; Bailin & Steinmetz 2005; Avila- 
Reese et al. 2005; Gottlober & Turchaninov 2006; Maccio et 
al. 2007; Bett et al. 2007; van den Bosch et al. 2002; Het- 
znecker & Burkert 2006; D'Onghia & Navarro 2007, e.g). 
Numerical simulations seem to agree with the early analyti- 
cal calculations, confirming that A's dependence on environ- 
ment is very weak (Faltenbacher et al. 2002). 

Problematically, the spin parameter cannot be mea- 
sured easily from observations, as none of J, E or M can 
be measured directly. As mentioned above, Jimenez et al. 
(2003) used an indirect method to determine the halo spin 
parameter, where they assumed a particular dark matter 
halo profile and an exponential disk for the baryonic matter 
and fit a parameterized model to observed rotation curves 
to recover A. Recently, Hernandez & Cervantes-Sodi (2006) 
(hereafter HC-S06) proposed another indirect method for 
approximating A, in the same vein as Jimenez et al. (2003) 
but without needing to use rotation curve information. 

Like Jimenez et al. (2003), HC-S06 assume the specific 
angular momentum of the dark matter and the baryons are 
equal, but they determine the galactic rotation velocities 
with the Tully-Fisher relation. They show from compari- 
son with numerical simulations that their estimator gives 
an accurate measurement of A (Cervantes-Sodi et al. 2007). 
Hernandez et al. (2007) applied this to the SDSS survey to 
determine the spin parameters of a volume limited sample of 
8, 000 color-and concentration-selected SDSS spiral galax- 
ies, finding that the spin parameter distribution P(A) does 
indeed follow a log-normal distribution in very good agree- 
ment with that predicted by analytic or numerical calcula- 
tions. In a more recent study (Cervantes-Sodi et al. 2007), 
the same group computed the correlation of A as a function 
of mass and of over-density of galaxies in the SDSS, find- 
ing an anti-correlation between A and galaxy mass and no 
correlation between A and environment. 

In the present study, by taking advantage of new results 
and the recent MOPED analysis that provides detailed star 
formation histories for galaxies in the SDSS-DR3, we modify 
and improve the HC-S06 method to measure A for a sam- 
ple of 53, 000 galaxies. With MOPED we can use stellar 
masses derived directly from the spectra and thus avoid us- 
ing the Tully-Fisher law twice as is done in HC-S06. We use 
our sample to study how the dark matter spin parameter 
shapes the star formation history of galaxies. We also inves- 
tigate how the spin depends on environment as measured 
by both the galaxies' distances to catalogued galaxy clusters 
and a A-weighted marked correlation analysis for comparison 
with the Faltenbacher et al. (2002) numerical simulations. 
We find that active star forming galaxies today have a sig- 
nificantly broader spin distribution, reaching to significantly 
higher values of A, than galaxies that formed most of their 
stars at z > 2.5. We also show that spin shows little depen- 
dence on environment, thus confirming the analytical and 
numerical results from hierarchical torquing theory. Our pa- 
per is organized as follows: in §2 we describe the sample of 
galaxies used and in §3 the method for approximating A. We 
present our main results in §4 and our conclusions in §5. 



2 SAMPLE 

Our original sample contains ~ 3 x 10 J galaxies for which 
MOPED star formation histories have been computed (Pan- 
ter et al. 2007). The galaxies were selected from the SDSS- 
DR3 Main Galaxy Sample (Strauss et al. 2002; Abazajian et 
al. 2005), satisfying 15.0 < r < 17.77 and spanning the red- 
shifts 0.01 < z < 0.25. For each galaxy, MOPED computes 
the star formation fraction (SFF) in 11 roughly logarith- 
mically spaced time bins, the metallicity in each of the 11 
bins, and the estimated E(B — V) dust extinction (which 
includes both internal and Galactic dust). MOPED also es- 
timates the present stellar mass M* of each galaxy, along 
with the total mass of stars that has ever formed in that 
galaxy. More about the MOPED algorithm and sample can 
be found in Panter et al. (2007) . 

Of those galaxies with known star formation histories, 
our sample includes only disk-dominated galaxies (SDSS 
fractional de Vaucouleurs weight fracDeV < 0.5) which are 
viewed face-on (SDSS observed axis ratio expAB > 0.6). 
While this latter cut is rather generous, MOPED's dust es- 
timate empowers a good correction for internal extinction. 
Unless otherwise stated, the results presented below persist 
for more stringent cuts on both of these parameters, albeit 
with higher noise. This subset contains ~ 5.3 x 10 4 galax- 
ies. Where necessary, photometric parameters from DR6 are 
used (Adelman-McCarthy et al. 2007). 



3 METHOD 

Since none of the quantities in Eq. 1 is directly measurable, 
some estimation is required. The more accurate method to 
determine A would be to employ galaxy rotation curves along 
with detailed modeling of galactic potentials (Jimenez et 
al. 1997; Gnedin et al. 2007), but constructing very large 
samples in this way is prohibitively difficult. HC-S06 have 
derived a rough estimator for A in terms of SDSS measure- 
ments alone; it relies on the Tully-Fisher relation (Tully & 
Fisher 1977) to provide disk velocities where necessary. 

We will give a quick overview of the HC-S06 method 
and will note where we modify their method to obtain more 
robust As. The HC-S06 model contains a baryonic galactic 
disk of mass Md with an exponential profile 

where r is the radial coordinate and Rd the scale radius. 
This disk sits in a truncated, spherical, isothermal halo of 
mass Mh which gives it a circular velocity Vd- In this simpli- 
fied model, the halo dominates the quantities J, E, and M; 
the disk is an observable tracer of them but contributes neg- 
ligibly to their totals. Our data limitations prevent us from 
any more detailed modeling of the dark halo NFW profile, 
but this toy model should be a good approximation, as long 
as galaxies have flat rotation curves. 

To determine J, the halo specific angular momentum 
j h — J/M h is taken to be equal to that of the exponential 
disk jd = 2VdRd- This is a very good approximation for 
the protogalactic state, where baryonic and dark angular 
momenta were well mixed, and is often assumed to hold 
for collapsed galaxies, as long as their formation histories 
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Volume Limited (z<0.055) 



BTF; a = 4.0 
MOPED stors; p=0.2 
MOPED stars + gas 




Figure 1. /'(A) for a volume limited sample with z < 0.055 and 
M r < —19.0 (chosen for ease of comparison with Hernandez & 
Cervantes-Sodi (2006)) for several estimates of Mh (solid colored 
lines) with the Hernandez & Cervantes-Sodi (2006) best fit over- 
plotted (dashed line). 



were smooth (Dalcanton et al. 1997), but its robustness is 
limited if baryons lose a significant fraction of their angular 
momentum to dark matter during collapse. This question is 
still an open one. Early simulations (Navarro & Benz 1990) 
indicated that the baryons retain only a small fraction of 
their specific angular momentum through collapse, but ? 
found baryonic spins to relate closely to dark matter spins, 
albeit with a considerable angular misalignment between the 
two. Recent studies (Kaufmann et al. (2007); Mayer et al. 
(2008)) claim that the problem is fraught with numerical 
difficulties, and that with sufficient resolution, the baryon 
component retains 80 or 90 percent of its specific angular 
momentum. Furthermore, physical effects such as supernova 
feedback effects may play a role. In spite of this uncertainty, 
it has been common in the literature (e.g. Mo et al. (1998); 
Jimenez et al. (1997)) to assume that the specific angular 
momentum of the baryons is the same as that of the dark 
matter. Assuming the halo to be virialized, E is simply half 
the gravitational potential energy of the halo: E — V d Mh/2. 

If this spin estimator accurately probes the dark matter, 
it should yield the redefined spin parameter A' introduced 
by Bullock et al. (2001), as a measure of spin that is inde- 
pendent of density profile. The mapping between A' and A 
depends on halo shape, and can be difficult in crowded sys- 
tems with poorly defined edges. Any disparity between the 
two is well within the expected scatter of the spin estimator, 
so we will refer to what follows, roughly, simply as A. 

In HC-S06, Md is determined from a general baryonic 
Tully-Fisher relation (BTF), 



M d = AtfV? 

where the constants a 



(3) 



633M0(fcms 1 ) . Such an Atf implies a disk mass 
for the Milky Way (V d = 220km/ s) of M d = 1O U M . 

A critical issue is that one needs to estimate the mass 
of the halo in order to estimate A. The best way to do this 
is not clear, so the procedure we have adopted is to inves- 
tigate three separate plausible relations between disk mass 
and halo mass. We find that all three give very similar quan- 
titative results, and qualitatively the conclusions are identi- 
cal whichever method is adopted. 

HC-S06 assume F = M d /Mh is a constant, which they 
place at 1/25, and estimate Mh for their galaxies. Nominally, 
the BTF should probe all the baryonic mass in the disk, 
but since V d is measured only indirectly through another 
stellar mass TF relation, it introduces a level of abstraction, 
with stellar light giving rotation velocity which in turn gives 
baryonic mass. 

Since MOPED provides a robust estimate of the disk 
stellar mass M», we can circumvent the need to rely on the 
BTF scaling relation and instead use more information from 
the real, observed galaxies. For high mass galaxies, M* m 
M d , but since our sample includes a range of masses, we 
need to map M* to M d in a consistent manner, taking into 
account the increasing gas fraction towards lower M d . We 
explore two methods to perform this mapping. 

First, we assume the same constant baryon-to-halo frac- 
tion F, but calculate M d = M* +M g from the MOPED stel- 
lar mass M» and an estimate for the associated gas mass M g . 
Calura et al. (2007) used the MOPED galaxies and Kenni- 
cutt star formation rate scalings to estimate the mean gas 
fraction (M g /M*) as a function of M*. We use the results of 
their Fig. 2 to calculate M g . Note that Calura et al. (2007) 
found some dependence on the fraction of baryons lost as 
a function of mass (so F is not constant). However, over 
the mass range of the vast bulk of the galaxies considered 
here, the dependence is not large, ±10%, and we ignore the 
variation here. 

Second, following Gnedin et al. (2007), we make no gas 
correction to M* but assume that Mh = M*/F, where F* 
is allowed to vary with galaxy parameters. Gnedin et al. 
(2007) suggest that F* (which they call m d ) should scale 
with stellar surface density as 



F, = F* 



m,r: 



10 9 - 2 M o fepc- 2 



(4) 



3.5 and Atf 



for some power p. As they observe, such a dependence on sur- 
face density follows naturally from a Kennicutt-like density- 
driven star formation law. By comparing observed and mod- 
eled TF and Fundamental Plane parameters and residuals, 
Gnedin et al. (2007) find a best ut p = 0.2. 

To summarise, to test our assumptions, we examine 
three assumptions for Mh'- (1) HC-S06's baryonic Tully- 
Fisher estimate ("BTF") where M h « A T FV d /F, (2) us- 
ing the MOPED stellar mass with the Gnedin et al. 
(2007) variable disk fraction ("MOPED stars") where M h « 
M»/F,(M,), and (3) using the gas-corrected MOPED stellar 
mass ("MOPED stars + gas") where M h » (M» + M g )/F. 
In all cases, we choose the constants such that the Milky 
Way has M d = IO^Mq and M h = 2.5 x 10 12 M Q . For the 
first case, we use a = 4.0 in close agreement with McGaugh 
et al. (2004), although the results are similar to the HC-S06 
choice of a — 3.5. 
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Finally, for each of these cases, we have an approxima- 
tion for the spin parameter A which looks like 

2 1/2 VjR d 
GM h ' 



BTF; a = 4.0 



A 



(5) 



Here, Rd is taken from the SDSS i-band expRad 1 , and the 
disk velocity Vd is determined by the i-band Tully-Fisher re- 
lation that log T/ d = -0.130(M 4 + 21.327) + 2.212 published 
by Pizagno et al. (2007) 2 . Mh is given by one of the three 
cases above. We choose i-band for Rd and Vd to limit po- 
tential confusion from active bursts of star formation, and 
use the MOPED dust parameter to correct for internal and 
Galactic extinction. Hernandez et al. (2007) used an older 
i?-band TF relation, which had to be converted to Sloan 
photometry. 

This estimator of A is very rough, and should not be 
trusted as precise for any one real galaxy. Typical errors, 
dominated mostly by scatter in the Tully-Fisher relation- 
ships, are of order at least 30% (Hernandez et al. 2007). 
With large numbers of galaxies, we still hope to build reli- 
able A distributions by averaging over the noise. More per- 
niciously, however, unknown systematic offsets in the co- 
efficients severely limit any statements that can be made 
about absolute values A. In deference to these limitations, 
our conclusions are drawn only from relative comparisons, 
and should be robust as long as the scalings hold. 

Cervantes-Sodi et al. (2007) test their version of this 
estimator against previously published N-body simulations, 
showing that it works reasonably well. Another ideal test 
would be to calculate detailed spin parameters for a large 
number of galaxies for which high-resoluation rotation 
curves are available, following Jimenez et al. (2003); Gnedin 
et al. (2007), comparing the results to the estimate based 
on SDSS parameters. Catinella et al. (2006) recently pub- 
lished a large sample of such rotation curves which that 
overlap significantly with the SDSS, but they contain very 
few galaxies at z > 0.055, where the MOPED analysis can be 
performed. The cumulative distribution of spins below do 
agree well with the results of Jimenez et al. (2003), but no 
truly solid test of these spin estimations has been performed 
so they must be interpreted with appropriate caution. 



M, > 4.9E+10 
Mtotol > 9.4E+10 
4.9E+10 > M, > 2.2E+10 
9.4E+10 > Mtotgi > 4.1E+10 
2.2E+10 > M, > 8.4E + 09 
4.1E+10 > M, oio i > 1.6E+10 
M, < 8.4E + 09 
Miotoi < 1.6E+10 




MOPED stars; p = 0.2 




MOPED stars + gas 




4 RESULTS 

First, we have tested that we produce a reasonable distri- 
bution of A's. Dark matter halos have long been predicted 
to exhibit a log-normal -P(A), with median A ranging from 
~ 0.03 to ~ 0.05 and logarithmic width <7io g A ~ 0.5 (Heav- 
ens & Peacock 1988; Vitvitska et al. 2002; van den Bosch et 
al. 2002; Hetznecker & Burkert 2006). 

Our P(A) is plotted in Figure 1 for the several Mh es- 
timates discussed above. For all assumptions, the distribu- 
tions are roughly lognormal and in general agreement with 
theoretical predictions and the HC-S06 best fit, which is 
overplotted. Remember that the coefficients going into Eq. 



1 In its exponential fitting, SDSS reports half-light radii R e which 
relate to our desired radii by Rd = il e /1.68. 

2 They actually publish Vgo, the velocity at a radius encompass- 
ing 80% of the i-band flux, which ought to be a good estimator 
of V d . 



Figure 2. For each estimate of M h (top - BTF, middle - MOPED 
stars, bottom - MOPED stars + gas), P(\) distributions for the 
entire sample split into quartiles by the present mass of stars 
in the galaxy M* (solid lines) and the total mass of stars ever 
formed in the galaxy M tota i (dotted lines). The curves are nearly 
identical, as would be expected if M, rises monatonically with 
Mtotai i an d the respective mass cuts are shown in the legend in 
the top panel. 



5 were independently set to reproduce the mass properties 
of the Milky Way, and although unknown biases in the es- 
timator could shift these distributions in either direction, 
they agree quite well. 

Next, we investigate how halo spin depends on MOPED 
stellar mass and star formation history. Fig. 2 shows P(A) for 
different stellar masses. Note that the distributions become 
narrower and with smaller mean as the stellar mass of the 
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galaxy increases. For stellar masses smaller than 1.5 x 10 
Mq the distribution extends to values larger than A = 0.1, 
while for higher masses the values of A end at A ~ 0.05. This 
trend is seen for all our estimates. In the MOPED stars 
case, increasing the surface density exponent p effectively 
decreases A for low mass galaxies, shifting their distribution 
to the left, and could potentially eliminate this mass depen- 
dence. The observed trend, however, persists for p < 0.5, 
which seems to be safely far above the Gnedin et al. (2007) 
best fit p — 0.2. This observed anticorrelation between mass 
and A was also recently observed by Cervantes-Sodi et al. 
(2007). 

Taking advantage of the previous MOPED analysis, we 
show in Fig. 3 how the spin distribution changes as a func- 
tion of the star formation history. The first panel compares 
-P(A) for galaxies which have formed more than 50% of their 
stars at look-back times older than 10 Gyr to that for galax- 
ies which have formed more than 1% of their stars in their 
last 0.2 Gyr. The next panel makes more extreme cuts for 
the star formation fraction in the same time bins (more than 
95% before 10 Gyr and 3% in the last 0.2 Gyr), and the next 
even more extreme (98% and 6% respectively). For com- 
pleteness, we show the three Mu estimates in each panel 
of Fig. 3. Among both the MOPED stars and MOPED stars 
+ gas curves, a very clear split is seen between narrowly 
peaked low-A, old stellar population galaxies and broadly 
distributed high-A, young stellar population galaxies. Fur- 
thermore, as the definitions of these regimes are made more 
extreme (higher SFF, lower panels), the differences between 
the A populations widen. The BTF curves are slightly more 
enigmatic, with the peaks shifting in the opposite direction, 
but with the MOPED trend persisting out in the tails. Since 
the MOPED stars and MOPED stars + gas methods more di- 
rectly probe the true mass of each galaxy we suggest those 
as representative of the underlying spin distributions. 

The recent star formation populations in Fig. 3 do ap- 
pear to roughly overlap with the low-mass quartiles in Fig. 2, 
in accordance with the finding by Heavens et al. (2004); Pan- 
ter et al. (2007) that galaxies with total stellar masses below 
5 x 10 10 Mq dominate star formation today. If star-forming 
galaxies are typically low mass, it appears from Fig 3 that 
they are also generally high spin. But furthermore, 30 — 60% 
of the galaxies contributing to the three recent star forma- 
tion curves in Fig. 3 have total stellar masses above 8 x 10 9 
Mq, suggesting that star formation is occurring today in 
high spin galaxies of all masses. 

Fig. 4 tries to pinpoint whether the correlation between 
spin and star formation seen in Fig. 3 is due primarily to 
link between mass and star formation or if spin provides 
a second parameter for star formation history. We define 
the "Recent SFF" as the fraction of stars formed within 
the past 0.2 Gyr, and plot it against stellar mass M* in 
Fig. 4. As discussed extensively in Panter et al. (2007), it 
shows that recent star formation occurs predominantly in 
low-mass galaxies. To determine if spin plays an additional 
observable role in galactic star formation, in the right panel 
we plot the residuals from the binned mean of the left plot 
as a function of spin for each of the three estimates. The 
best fit lines for the residuals, of the form 



SFF before 1 Gyr > 0.50 
SFF since 0.2 Gyr > 0.01 




SFF before 1 Gyr > 0.95 
SFF since 0.2 Gyr > 0.03 





log(Recent SFF) - mean(M,) = alog(A) + b, 



(6) 



-1.1,6 



-0.8 with a scatter of <7R CCCnt sff = 0.47 



Figure 3. For each estimate of (dotted lines - BTF, dashed 
lines - MOPED stars, solid lines - MOPED stars + gas), P(X) is 
shown for subsamples split on the basis of the fraction of stars 
formed before 10 Gyr in lookback time (red lines) and the fraction 
of stars formed since 0.2 Gyr in lookback time (blue lines). The 
star formation fraction (SFF) cuts are shown in each panel, and 
lower panels show greater extremes of the population. From top 
to bottom, each curve is calculated from roughly 15, 000, 2, 500, 
and 1 , 000 galaxies respectively. 



(BTF), a = 0.92, b = 0.61 with cT Rccell t sff = 0.49 (MOPED 
stars), and a — 0.50,6 = 0.29 with o"R CCGn t sff = 0.50 
(MOPED stars + gas). That the BTF residual (top) shows a 
negative slope is consistent with the enigmatic result seen in 
Fig. 3, where by any measure the correlation between A and 
SFF is less clear. Both the MOPED stars (middle) and MOPED 
stars + gas (bottom) show slight positive slopes, indicat- 
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ing that high-A galaxies for a given M» host an excess of 
recent star formation. Although the scatter is considerable, 
Pearson's p has a virtually zero null probability in each case. 

The discrepancy between the three A estimates is trou- 
bling. Since we expect that the MOPED masses more accu- 
rately reflect true galaxy masses, as they are not plagued by 
the BTF's double use of the same scaling relation, Fig. 4 could 
indicate a slight dependence on spin as a second parameter 
governing the star formation history of galaxies. However, 
because ultimately we sample only the luminous disk matter 
whose link to halo properties may be very complicated, the 
question of spin an independent drivers of star formation 
history remains poorly resolved. 

4.1 Enviromental dependence 

We will now investigate the role of environment in deter- 
mining spin. With the largest sample of (indirectly mea- 
sured) spins yet collected, we are well positioned to answer 
this question. We approach this question in two ways: first 
we compare to a sample of galaxy clusters found in the 
SDSS-DR6 catalogue using a new algorithm (Berta 2008, in 
preparation), and second we use a marked correlation study 
(Sheth et al. 2006) to compute the clustering properties of 
the spin itself. 

The cluster catalogue was recently constructed from 
DR6 galaxy photometry with the Rosy Matched Fil- 
ter (RMF) cluster-finding algorithm, which integrates a 
Matched Filter technique (Postman et al. 1996; Kepner et 
al. 1999; Kim et al. 2002) with an early-type red-sequence 
search criterion like that used by the Red-Sequence Cluster 
Survey and the maxBCG cluster catalogue (Gladders & Yee 
2000; Koester et al. 2007). Clusters are identified as peaks in 
a 3D likelihood map measuring how well observed galaxies 
match a realistic cluster filter in a 5D position-brightness- 
color-color-color space, and characterized by a richness pa- 
rameter Armf, designed to scale with total cluster lumi- 
nosity. Note, the RMF catalogue was compiled only from 
bulge-dominated galaxies with f racDev > . 5, a completely 
disparate set from the galaxies studied here. 

For each galaxy for which we measure a spin, we cal- 
culate a radius Rcyiinder such that an imaginary cylinder of 
depth Az = 0.02 and radius R cy iinder is just large enough 
to include the one RMF cluster. This is basically the pro- 
jected distance to the nearest cluster, but it is expresed in 
this manner because precise galaxy-to-cluster distances are 
nearly impossible to measure (Finger-of-God effects, cluster 
redshift errors, etc.). We will use Rcyiinder as a rough tracer 
of the local density near each galaxy, with smaller Rcyiinder 
corresponding to denser environments. Fig. 5 shows the cu- 
mulative distribution of Rcyiinder for high and low spins, 
with each curve being normalized to its last bin. In these 
plots, shallower slopes correspond to more clustered galax- 
ies. 

The first panel of Fig. 5 seems to indicate that low-A 
galaxies are much more strongly clustered than high-A galax- 
ies, and such a finding is in accordance with Avila-Reese et 
al. (2005), but is this difference a result of the spins alone? 
It has long been known that clusters are overly populated by 
massive, early-type galaxies with old stellar populations. As 
indicated above, these galaxies are more likely to sit in low-A 
halos and could thus contribute to the observed clustering 
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Figure 6. A marked correlation function showing the A-weighted 
galaxy pair counts as a function of separation radius WW(r) 
normalized by the unweighted pair counts DD(r). Three volume 
limited samples are shown (z < 0.055, z < 0.1, z < 0.15). Error 
bars are calculated from 6 jacknife realizations of WW(r)/ DD(r). 



difference. The second panel shows only low-mass galaxies 
{Mtotai < 5 x 10 10 M , roughly the median of the com- 
plete sample), and as expected, the difference in clustering 
strengths decreases, although still indicating low-A galaxies 
to lie preferentially nearer to clusters. In the third panel, we 
include only the low-mass (M tota i < 5 x 10 10 Mq) galaxies 
which have formed more than 1% of their stars since 0.2 Gyr 
(this corresponds to top panel of Fig. 3), and all visible clus- 
tering dependence on spin vanishes. It seems the observed 
Rcyiinder — A correlation can be adequately explained as the 
confluence of galaxy mass and stellar age increasing towards 
cluster centers. 

Using RMF clusters as density tracers has the advan- 
tage of covering a wide volume and broad range of galaxy 
masses, but to investigate further the correlation proper- 
ties of the dark halo spin we performed a marked correla- 
tion analysis (Sheth et al. 2006) in several volume limited 
samples, using galaxy spin A as the mark. Fig. 6 shows the 
deviation of A from its mean value as a function of scale, 
as calculated in three volume limited samples: roughly the 
HC-S06 volume z < 0.055, roughly the "faint" volume from 
Sheth et al. (2006) z < 0.1, and a deeper volume z < 0.15, 
using the MOPED stars + gas estimate in all cases. The as- 
sociated absolute magnitude cuts are shown in Fig. 6. 

We find no significant scale dependence of A, a re- 
sult which persists among the BTF, MOPED stars, and MOPED 
stars + gas estimates. Each volume corresponds to a much 
narrower range of galaxy masses than Fig. 5, which accounts 
for the absence of any such clustering dependence as seen in 
the first panel of Fig.!5. 

This is what is predicted in analytical and numerical re- 
sults (Heavens & Peacock 1988; Nagashima & Gouda 1998; 
Bett et al. 2007) but in contrast with the results of Fal- 
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Figure 4. Le/t: Scatter plot showing the recent star formation (< 0.2 Gyr) as a function of stellar mass M*, with the best linear fit 
(dashed line, both panels) and the mean in bins of AlogM* = 0.1 (solid line, both panels). Right: The residual, having subtracted the 
binned mean (solid line) from each point at left, as a function of spin for three halo mass estimates: BTF (top), MOPED stars (middle), 
MOPED stars + gas (bottom). 
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Figure 5. In each panel, the fraction of galaxies for which at least one RMF galaxy cluster is found within a cylinder of depth Az = 0.02 
as a function of the radius of that cylinder R C yiinder- High spin (A > 0.04, blue lines) and low spin (A < 0.04, red lines) are plotted for 
the entire sample (left), low-mass galaxies (center), and low-mass, star-forming galaxies (right). Poisson errors are shown as lower limits 
on the true error. Note that these plots are cumulative; the individual points are not independent. 



tenbacher et al. (2002) who explicitly computed the marked 
correlation of the spin parameter A in their simulations and 
found an environmental dependence. However, Faltenbacher 
et al. (2002) argue that large samples are needed to measure 
this signal, it therefore remains to be seen if much larger 
sample can confirm our result. Cervantes-Sodi et al. (2007) 
also recently found no environmental dependence for A. 



5 CONCLUSIONS 

We have used an indirect method to compute the dark mat- 
ter spin of galaxies and combined it with MOPED determi- 
nations of the star formation history of the SDSS galaxies. 
Exploiting the large sample available, we have studied the 
influence of dark matter spin on galaxy mass, star formation 
history, and environment. We find that galaxy dark matter 
spin and stellar mass are anti-correlated: lower stellar mass 
galaxies exhibit broader and generally higher distribution 
of spins than high-mass galaxies. Furthermore, according 
to halo mass estimates determined from the MOPED stellar 
masses, galaxies which have formed 1 — 6% of their stellar 



mass in the past 0.2 Gyr also have typically broader and 
higher-A spin distributions than galaxies that have formed a 
large fraction of their stellar mass at look-back times larger 
than 10 Gyr. Although mass is the prime parameter de- 
termining the current star formation rate, the galaxy spin 
parameter might play a weak secondary role, with higher- 
spin galaxies having more current star formation at given 
stellar mass. 

Such a result could be explained in terms of high-A ha- 
los hosting low surface density disks with consequently less 
efficient star formation rates in the past leaving a larger 
reservoir of gas today. This role of spin in galaxy formation 
was also emphasized by Jimenez et al. (1997), who proposed 
a 'spin bias'. We find here that though this effect might exist 
it is likely a subdominant effect. Along the lines of Vitvit- 
ska et al. (2002); Mailer et al. (2002), the link between star 
formation and A might follow from high spin galaxies expe- 
riencing recent mergers. Even if these mergers are merely 
the signature of the large scale tidal torques on hierarchical 
galaxy formation, they might both increase spin and spur 
starbursts in an individual galaxy, provided the stellar disk 
survives the merger. 
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We have also looked at environmental effects: using the 
RMF catalogue of galaxy clusters in the SDSS we find a 
very weak anti-correlation between the value of dark matter 
spin and proximity to a cluster, but such as would be con- 
sistent with mass and star formation being positively cor- 
related with cluster proximity (see Mateus et al. (2008) for 
a comparison of the environmental dependence of galaxies 
with numerical models of galaxy formation). A marked cor- 
relation study also shows no strong correlation with galaxy 
separation, in agreement with previous studies' predictions. 

Finally, we note that this analysis requires an estima- 
tion of three unobservable galaxy parameters M, J, and E, 
which is not at all a trivial problem. Determining the latter 
two can be difficult even in numerical simulations: although 
each dark matter particle can be followed, arbitrarily de- 
fined galaxy borders can result in arbitrary results for A 
(D'Onghia & Navarro 2007). Furthermore, as the link be- 
tween the observed stellar disk and the total halo angular 
momentum is complicated (van den Bosch et al. 2002, ie), 
the problem is very difficult. We used three estimates for 
the halo mass Mh ~ M, and found that they did change our 
results somewhat. 

Although MOPED probably gives a much more accurately 
measure of the true baryonic mass of a galaxy than the 
double-scaling relation BTF, further study is necessary to 
elucidate the accuracy of the estimator used here. 
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